Localization and circulating currents in curved graphene devices 
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We calculate the energy spectrum and eigenstates of a graphene sheet which contains a circular 
deformation. Using time-independent perturbation theory with the ratio of the height and width of 
the deformation as the small parameter, we find that due to the curvature the wavefunctions for the 
various states acquire unique angular asymmetry. We demonstrate that the pseudo-magnetic fields 
induced by the curvature result in circulating probability currents. These circulating currents in 
turn produce local real magnetic fields ~ 100 /iT which can be measured using current technology. 



Graphene, a monolayer of carbon atoms in a honey- 
comb lattice, has attracted a lot of interest in the last 
decade [TH3]- Along with many interesting electronic 
transport properties [Ij graphene has some intriguing 
structural and mechanical properties [S]. A suspended 
graphene sheet can be deflected by applying a gate volt- 
age [6j. Early on, curvature induced by topological de- 
fects was studied in the context of carbon nanotubes 
and fullerenes [7HTD]. Graphene sheets deposited on a 
substrate naturally show corrugations [TTJ [T^ which can 
be modeled as defects [13]. Deformations in a graphene 
sheet can provide a mechanism for spin relaxation which 
is stronger than relaxation due to intrinsic spin orbit in- 
teractions [Uj. These deformations can also lead to sig- 
nificant effective pseudo-magnetic fields [TS] . Two recent 
works have investigated properties of pseudo-magnetic 
fields in graphene sheets: In Ref. [Hj de Juan et al. in- 
vestigate Aharonov-Bohm interferences in the local den- 
sity of states due to a fictitious strain-induced magnetic 
field, and in Ref. [T7] Abedpour et al. calculate the 
pseudo-magnetic field induced by shear stress in circular 
graphene rings. 

The purpose of this Letter is to investigate the possi- 
bility of exploiting curvature, induced by elastic defor- 
mations, in a graphene sheet for generating states with 
circulating currents. We focus on systems where the cur- 
vature in the graphene sample can be controlled exter- 
nally. A possible way of realizing such systems is to sus- 
pend graphene samples on a substrate where an annu- 
lus has been etched away. Modeling the curvature as a 
Gaussian-shaped bump and using perturbation theory, 
we show that it is possible to generate localized circulat- 
ing carriers in the ground state for such systems. 

According to the standard theory of elasticity [131 HBl" 
[5D] the effect of curvature can be captured by including 
a scalar (Vi) and vector (V2 — — i-^ps) potential 
in the Hamiltonian jl3j with Vi = gi(uxx + Wyy) and 
V2 = g2{uxx - Uyy + 2iuxy). Here Uij (with i,j G {x, y}) 
is the usual deformation tensor, 
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FIG. 1. (Color online) (a) Schematic of type (a) deformation, 
(b) Schematic of type (b) deformation, (c) Pseudo-magnetic 
field produced by type (a) deformation, (d) Pseudo-magnetic 
field produced by type (b) deformation. See the text for fur- 
ther details. 

with Ui = Ui(r), i S {x,y,z}, the Cartesian components 
of the average displacements. The coupling constants 
gi and 52 are known from transport measurements |13j . 
We assume that the length scales associated with the 
deformation are much larger than Occ, where Occ ^ f.4 
Arepresents the lattice constant of graphene, so that in- 
ter valley scattering can be neglected. The single-valley 
Dirac Hamiltonian for the system in the absence of ex- 
ternal magnetic fields is given by: 

-Tj f -eVi px - ipy + eV2 \ 

\Px + Wy + eV^ -eVi J 

Here vp is the Fermi velocity and p—~ih\/ represents 
the momentum operator. Now we focus on the shape of 
the deformation. In particular we consider two deforma- 
tions: (a) a Gaussian centered at the origin, i.e ^^(r) = 

Ae and (b) a ring-shaped Gaussian deformation 

modeled as Uz[r) — Ae = 

(Vx2+„2_,)2 

Ae . These deformations are schematically 

depicted in Fig. 1(a) and 1(b). We consider deforma- 
tions with small height and large width and apply per- 
turbation theory, using a = A/b as the small parameter. 
In polar coordinates the potentials Vi and V2 are given 
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FIG. 2. (Color online) Contour plot of ^VPfe^ , (a) for 
the lowest electronic state of type (a) deformation, (b) for the 
first excited state of type (a) deformation, (c) for the lowest 
electronic state of type (b) deformation and (d) for the first 
excited state of type (b) deformation. Parameters used are 
a = 0.04, 13 = 0.1, 7 = 0.4, egi = 3 eV, eg2 = 2.276 eV [l3]. 
See the text for further details. 



by l/i = 2gij^fir) and V2 = 2g2=^e*'V(r-), where 
Tn-e^^^^/''^ for type (a) deformations 
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The Hamiltonian [Eq. (2)] can be written in polar coor- 
dinates as "H = "Ho + Hi with 



and 
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Here Hq describes the unperturbed system, which is 
a fiat circular graphene sheet of radius R and Hi is the 
perturbation. The normalized eigenstates of the unper- 
turbed Hamiltonian Ho are given by : 



^'fc,m(r,6i) = 
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where Jm(rfc) is the Bessel function of integer order m 
with energies em,k — shvpk. s = ±1 represents the con- 
duction and the valence band. It should be noted that 
the states described by Eq. (6) are also eigenfunctions of 
the z-component of the total angular momentum L^'- 

L^-^k,ra = i-ihd^ + ^<T3)*fc,m = (m + (7) 

Now we impose infinite mass boundary conditions |22| 
which are valid in the regime R 3> acc- This can be 



done individually for each m and results in the eigenvalue 
equation: 



Jmikm.lR) — Jni+l{km,lR) : 
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where the integer I represents the ascending order of the 
values I k„ij \ that satisfy the above equation. The eigen- 
values for m and —to — 1, which correspond to equal and 
opposite total angular momenta, are not degenerate. The 
spectrum resulting from Eq. (8) is discrete due to the 
finiteness of the graphene sheet, and exhibits a band gap 
at the Dirac point {k^.i = 0). The solution of Eq. (8), 
closest to the Dirac point is given by fc_i__i = — ^'^^^"^ 
and fco.i = -^^^^^ The magnitude of the band gap is thus 
2nvF^^^ ~ 30 meV for i? = 50 nm ;23;. It is useful to 
write the perturbation Hamiltonian "Hi as 



17 2ea2 



526 -51 



(9) 



where fi = b/R and x — r/R. The form of the deforma- 
tion is now given as 



x^e for type (a) deformations 



[x — 7)^6 



^ for type (b) deformations. 

(10) 

Here 7 = c/i? is a measure of the location of the de- 
formation for type (b) deformations. The vector poten- 
tial in our formalism generates a pseudo-magnetic field 
Bps = V X Aps, perpendicular to the graphene sheet. 
For type (a) deformations this magnetic field is Bps'* = 
852^46" 6^ sin(36')z as sketched in Fig. 1(c). Similarly, 



for type (b) deformations Bp, 
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2r — 2c r — 6cr^ + b e — b cr) sin(30)z as shown in Fig. 
1(d). The pseudo- magnetic fields are not real magnetic 
fields but a consequence of the deformation. Although 
these fields are zero when averaged over the entire circu- 
lar graphene sheet (as can also be seen from Figs. 1(c) 
and 1(d)), they can locally be quite large, see the discus- 
sion at the end of the Letter. 

Now we calculate the corrections due to the curvature 
to the energy spectrum and eigenfunctions up to second 
order in a^. The energy eigenvalues are then given by 



.(0) 



are given by 



(1) 



Similarly, the wavefunctions 
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eigenvalues e^. 
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The unperturbed 
and wavefunctions ^?'' are given by 



JO) 



shvpkm I and 



(0) 



(11) 



isJ™+i(rfc™j)e*^"+')'' 
where Ami is the normalization factor given by 

Ami = 2 / [Jm{km,ix) + Jm+l{km,ix)] xdx. (12) 
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The corrections e[.^ ^ , e^'^ ^ and ^'[.^ ^ are calculated by 
standard time-independent non-degenerate perturbation 
theory [53] and we obtain: 



(1) 4ea^gi f{x) [j^^^jk^jx) + J?n+i{km,ix)] xdx 
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Here the matrix element Vk^, i,;km i is given by: 
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The angular integral in Eq. ( 14 ) is only non-zero if m = 



m , m = m 
calculations. 



3 or 



3, which simplifies the 



Now we present our numerical results. In principle the 



and 4-^^^ [Eqns. (|l3l)] 



sum in the expressions for el^^ 

runs over an infinite number of values of km,h but we 
restrict the sum to values of k„ij such that m £ (—6, 5) 
and I S (—15, 15). We have verified that the results ob- 
tained by using this restricted sum presented for the low- 
est two electron-like and lowest two hole-like eigenvalues 
and eigenfunctions are converged up to at least four sig- 
nificant digits. 

In Fig. 2 we plot the probability density ^j^^ 
over the entire graphene sheet for the lowest electronic 
state (i.e km.i = 1.4347/i?) and the first excited electronic 
state (i.e km,i = 2.62987/i?) for both type (a) and (b) 
deformations. A contour plot of i is shown in 

Fig. 2(a) for the lowest electronic state (ground state) for 
a central Gaussian [type (a)] deformation. We see that 
there is no angular dependence of the probability density, 
which is in agreement with previously obtained results 
[mill]. In contrast. Fig. 2(b) shows that ^^-fc^,, for 
the first excited electronic state of a type (a) deformation 
is angle-dependent. This angular dependence is in agree- 
ment with the pseudo-magnetic field B^"^ generated by 
such a deformation. For a ring-shaped Gaussian, on the 
other hand, it is evident from Figs. 2(c) and 2(d) that 
the probability density ^'^^ is angle-dependent for 

both the ground state and the first excited state. The 
fact that also the probability density of the ground state 
is angle-dependent makes the ring-shaped deformations 
more suitable for experimental demonstration. It should 
be noted that the probability density is high around the 
location of the deformation (r = OAR). As we increase 
a, the effect of deformation-induced localization of the 
wavefunctions for both types of deformations increases. 
A localized probability density results in circulating cur- 
rents which in turn generate local magnetic fields. In 




FIG. 3. (Color online) (a) Vector field plot of j over the en- 
tire graphene sheet for the lowest electronic state for a ring- 
shaped Gaussian [type (b)] deformation, (b) Classical trajec- 
tory in real space calculated using Eq. (15). Solid (black) 
and dashed (red) lines represent two different sublattices. (c) 
Classical trajectory in momentum space calculated using Eq. 
(15). Again solid (black) and dashed (red) lines represent two 
different sublattices. See the text for further details. 



order to investigate this effect we calculate the probabil- 
ity current density j = f_F^'^ ,'''^fc„,- Fig. 3(a) shows 
a vector field plot of j over the entire graphene sheet 
for the lowest electronic state of a ring-shaped Gaussian 
[type (b)] deformation. 

Localization of the carriers is a repercussion of the 
topological effect of the pseudo-magnetic field. To eluci- 
date this further we plot in Fig.[3]the classical trajectories 
calculated from the corresponding classical Hamiltonian 
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Fig 3(b) shows the position space cyclotron orbits for 
the ground state of the ring shaped Gaussian defor- 
mation [25]. The two different sublattices of graphene 
have different orbits as is shown by the sohd and the 
dashed hues. Fig. [sj^c) shows the corresponding momen- 
tum space closed orbits. These classical orbits exhibit 
remarkable similarity to the quantum-mechanical proba- 
bility density and the probability density current shown 
in Figs. [2]Jc) andjsjja). These circulating orbits produce 
local magnetic fields (which are zero when averaged over 
the entire graphene sheet). It should be stressed that the 
existence of these magnetic fields does require the carri- 
ers to be valley-polarized, since otherwise the degenerate 
carriers in the other valley will cancel these local mag- 
netic fields exactly. As has been suggested recently, a 
special kind of line defect can be used in graphene as a 
valley filter [26j. Alternately, weak intervalley scattering 
has been shown to lift the valley degeneracy [16]. Ex- 
perimentally one can measure the angular distribution of 
the probability density using a scanning tunneling micro- 
scope (STM) . To get an estimate for the order of magni- 
tude of the magnetic fields produced we approximate the 
current density depicted in Fig. [sjja) by a circular wire 
of radius OAR, where six equal arcs on the circumference 
carry alternating counter circulating constant currents. 
For typical carrier densities of 10^^ cm~^ the resulting 
magnetic fields are then of order 100/iT. Magnetic fields 
of this order can be easily measured by using a single 
nitrogen- vacancy (NV) impurity in diamond [27j . 

To summarize, we have calculated the energy spectrum 
and wavefunctions for two circular Gaussian shaped de- 
formations in a graphene sheet using perturbation theory. 
Due to the curvature in the graphene sheet the carri- 
ers become localized around the deformation and follow 
circulating orbits. These circulating carriers produce lo- 
cal magnetic fields which are zero when averaged over 
the entire sample. Interesting angular dependence of the 
probability density is predicted for the ground state of 
a ring shaped Gaussian deformation. We suggest differ- 
ent ways to experimentally verify our predictions with 
current technologies, which will provide further insight. 
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